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Resonant scattering of muonic hydrogen atoms via back decay of molecular com- 
plex, a key process in the understanding of epithermal muonic molecular formation, 
is analyzed. The limitations of the effective rate approximation are discussed and 
the importance of the explicit treatment of the back decay is stressed. An expression 
of the energy distribution for the back-decayed atoms is given. 
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Direct measurements of epithermal resonant molecular formation, recently 
reported by the TRIUMF Muonic Hydrogen Collaboration required de- 

tailed considerations of processes which were not previously well-studied. Reso- 
nant scattering of muonic atoms via back decay of the muonic molecular complex 
(MMC) is one such process. Despite substantial theoretical efforts in improving 
the accuracies of muonic hydrogen scattering cross sections, so far little attention 
has been paid to the back decay process as a scattering mechanism of the muonic 
atom (except for the spin flip in the dfid system). It is the purpose of this paper 
to emphasize the importance of the resonant scattering and the associated MMC 
dynamics with the hope of stimulating further theoretical studies. 

Resonant formation of MMC, ^xa^ +DXy-Ki — > [(d^a)f^xee],yj,i<-j,taking place 
with the rate ^^^k u K i-^t^a), is generally followed by competing processes of 
either stabilization with the effective rate A/ leading to fusion, or back decay 
to /ia^ + DX^'x' with the width ,r^,. Here a = d,t, x = p,d,t, and 

X = H, D,T. F is the hyperfine state of fia, S is the spin of dfj,a, and ViKi, 
i^fKf are the vibrational and rotational quantum number of DX and MMC 
respectively. Our notations follow those of Ref. but we explicitly account 
for the possibility of (de)excitation of DX upon back decay i^lK- ^ I'iKi, i.e., 
resonant (de)excitation, whose importance will become apparent below. 

In the analysis of conventional /xCF experiments, an effective renormalized 
formation rate has been widely used, into which the effect of the back decay 
is absorbed as: 

>'diia= "^ufKfY^K.^ljK^vjK;^ (1) 

where W^K^ = Xf/(\f + Y.u'k' ^ffxt v'K') fusion probability, and uok^ is 

the initial Ki population. We observe that even in the case where the transport 
of muonic atoms can be neglected, at least one of the following criteria must met 
in order to justify the effective rate approximation of Eq. |l| in describing fusion 
yields: (a) trivial condition that the back decay probability (1 — W^^^^) <C 1, (b) 
rapid (compared to MMC formation) re-thermalization of fia in the equilibrium 
condition, or (c) negligible change in fia energy in lab frame before and after the 
back decay. For example, the condition (a) is satisfied for dfit formation at low 
energies, while the condition (b) applies for dfid (at least at high densities), and 
(c) may be possible in condensed matter if recoil-less processes dominate. If none 
of the above criteria are satisfied, the fj,a can be removed from the resonance 
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regions affecting tlie kinetics, which was indeed the case in the experiment of 
Refs. mH,^. In simphfied model calculations, where one interaction, either of 
resonant scattering or potential scattering, is assumed to be sufficient to remove 
fia from the resonance region, the effective rate approximation overestimates 
the fusion yield by (A^^.^^^j^^ + Xscat)/iW^fKf^u5<i,ufKf + ^scat) where \scat 
is the potential scattering rate, compared to the explicit inclusion of resonance 
scattering channel [|| . For ij.t + D2 at resonance peak energy, this factor is as 
large as ~ 1.5. 

A more realistic estimate of the effect of resonant scattering requires the 
accurate /ia energy distribution after back decay, which depends on the details of 
the MMC dynamics including: (i) MMC recoil from fia impact upon its forma- 
tion, (ii) thermalization of MMC center of mass motion in collision with the tar- 
get medium, (iii) collisional relaxation/excitation of MMC ro-vibrational states 
{vfKf —>■ VjKj), (iv) MMC decay with possible DX excitation, and (v) DX 
recoil upon MMC decay. 

Apart from trival, but often a neglected effect of the kinematics (i,v), the 
cross section for [{diit)dee\ + D2 elastic collisions has been calculated by Padial et 
al. Q for relatively low energies; its extrapolation to epithermal energies suggests 
a value of about 3 x 10~^^ cm^. The corresponding collision rate ~ 2 x 10^^ s~^ is 
an order of magnitude larger than the MMC decay rate 10^^ s~^, hence substantial 
thermalization can be expected at high densities. 

Calculations of rotational transitions in the MMC were reported by Ostro- 
vskii and Ustimov |^] , and by Padial et al. |^ for the case of thermal equilibrium 
targets. Ostrovskii and Ustimov estimate relaxation rates of the order of 10^'^ 
s~^ (a value used in Ref. Q), while Padial et al., claiming higher accuracy, give 
~ 0.3 X 10^'^ s^^ at 300 K (rates are normalized to liquid hydrogen density). 

To date, there are no accurate calculations available for vibrational quench- 
ing of MMC, except for a rough estimate by Lane who gives 10''^ s~^ at room 
temperature, a rate much slower than other processes. He predicts, however, 
increasing rates for higher temperature and increasing v. Future accurate cal- 
culations of this process is highly desirable, as the consequence of non-negligible 
quenching would also have a significant impact on other aspects of molecular 
formation, such as the fusion rate A/. 
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Given the initial condition Vi,Ki,E^a, the approximate energy distribution 
of back decayed /io, f{E'^^), can be estimated from the following expression |^]: 



u'K'S v',K' 

X j dE^D{Eu)I [E'^a - El^K^VfKfS- u[K'^ - E^] , (2) 



where 



K.K, i^fKfS] v[K'i) = [eresivfKf) - AE^^K^y^K,^ (3) 

is the decay Q-value for the specific channel {vf,Kf,S) {v[,K[^S), with 
A.Ey.j^.y j^i = — [E{u[K'j) — E{viKi)\ being the binding energy difference be- 
tween the initial and final state of DX [i.e., (de)excitation energy of DX due to 
the resonant scattering). /[A] is the resonance intensity profile for detuning A 
{e.g., /[A] = (5(A) in the classical Vesman model), while D{Eu) is the Doppler 
broadening distribution due to the motion of MMC at the time of back decay 



{e.g., Gaussian distribution with the width a = ^4kTM^a/MMMC for thermal- 
ized MMC iiE'^^ > kT). 

9ufKj {^i, - (4) 

u',K' 

is the branching ratio for the decay into the state {v[,K[,S)., given the MMC 
state of {vf,Kf,S), and h^j^.{vfKfS\E^a) is the ro-vibrational population of 
MMC at the time of back decay, given the initial condition Ui,Ki and E^a- A 
full evaluation of Eq. ^ would require, in addition to the back decay matrix 
elements, a solution of kinetics equation involving all the competing processes; 
some limiting cases were considered in Ref. ||]. We note that in Refs. |10|,11|, 



back-decayed is suggested to have a thermal energy distribution of the target 
temperature, but we find this not be the case, even if MMC translational motion 
and/or ro-vibrational states are completely thermalized (see EqJ^. 

In our analysis for Ref. E'^f was varied between 1 meV to 0.3 eV in Monte 
Carlo calculations ||l2| explicitly taking into account the resonant scattering (an 
improved version of earlier calculation ||l^ ) to phenomelogically investigate its ef- 
fect, and some 7% difference in fusion yield was observed, giving a non- negligible 
contribution to the total systematic uncertainties. Although this error is not 
overwhelming, a substantial improvement in the accuracy of resonant formation 
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measurements would require, among others, detailed understanding of resonant 
scattering processes and MMC dynamics, the first step of which has been illus- 
trated in this report. 
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